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Abstract 

In the present work we use three-body interaction formalism to investigate the A'-multi-p interactions. 
First, we reproduced the resonances A(1270) and Ari(1270) in the pp and pK two-body interactions respec¬ 
tively, as the clusters of the fixed center approximation. Then, we study the three-body K — pp{f 2 ) and 
p — pK{Ki) interactions with the fixed center approximation of the Faddeev equations. Furthermore, we 
extrapolate the formalism to study the four-body, five-body and six-body systems, containing one K meson 
and multi-p mesons. In our research, without introducing any free parameters, we generate the A'2(1770) 
state in the three-body interaction with the mass of 1707 MeV and a width about 113 MeV, which are 
consistent with the experiments. We also find a clear resonant structure in our results of the five-body 
interaction, with a mass 2505 MeV and a width about 32 MeV or more, which is associated to the iF4(2500) 
state, where we obtain consistent results with the experimental findings. Furthermore, we predict some new 
states in the other many body interactions, 1 ^ 3 ( 2080 ), iF5(2670) (isospin / = 1/2), and Ar4(2640) (isospin 
I = 3/2), with uncertainties. 

PACS numbers: ll.80.Jy; 12.38.Lg; 13.75.Lb; 12.39.Fe 


1 


I. INTRODUCTION 


To understand the nature and structure of the particles found in the experiments, and to search 
for new hadronic states, are the main issues in nowadays particle physics which has been accepted 
generally that quarks are the basic building blocks of matter. Within the Gell-Mann-Zweig quark 
model [l|, for the normal hadron states, mesons are made of a quark-antiquark pair, qq, and 
baryons made of three quark components, qqq. On the other hand, there are some states found 
in the experiments, such as the mesons, /o(500), /o(980), ao(980), k(800), the baryons, A(1405), 
A^(1440), A^(1535), with structure and properties difficult to explain by the normal quark model. 
These states might therefore be called “exotic” states (recent experimental discussions are given 
in Refs. BB). To understand the structure and properties of these exotic states in the strong 
interaction, we need to exploit other theories or approaches, for example chiral perturbative theory 
Ba, effective held theory Lattice QCD QCD sum rule [isl- liil l , D yson-Schwinger 


equations 22H24i|. chiral quark model 25l427l|. chiral unitary approach (ChUA) [28l433l| . and so on. 
But, some discovered particles, such as the (?!)(2170) (also called X(2175), y(2175)), T(4260), 
A^*(1710), look like having a more complicated structure and could come from multi-body hadron 

physics drawing much attention for a long time 33-38|. 
1391 ] develops the ChUA for the three-body interaction, 


interaction, which is a subject in hadron 
With this motivation, the work of Ref. 
which combines the three-body Faddeev equations with on shell approximation of the ChUA and 
has reported several S'-wave resonances qualifying as two mesons-one baryon molecular 

states. In Ref. [i^ this combination of Faddeev equations and chiral dynamics in the DKK 
system obtains consistent results with QCD sum rules. When in some cases there are resonances 
(or bound states) appearing in the two-body subsystem of the three-body interaction, Ref. [Uj 


takes the fixed center approximation (FCA) [3a, I42l-l46l| to the Faddeev equations, where several 


multi-/9(770) states are dynamically produced, and the resonances /2(1270)(2'''+), /93(1690)(3“ 
/4(2050)(4++), P5(2350)(5 ), and /6(2510)(6'''+) are theoretically found as basically molecules of 

an increasing number of /9(770) particles with parallel spins. Analogously, the resonances ATKlddO), 
iFg (1780), Ar|(2045), 775(2380) and a new ATg are produced in the i7*-multi-/? systems and could be 
explained as molecules with the coniponents of an increasing number of /?(770) and one AT* (892) 


47t |. Also in Ref. charmed resonances ZI 3 , D 


4 ! 


Dg and are predicted 


meson in Ref. _ _ ^ 

in the ZI*-multi-/? interaction. Note that, recently a resonance structure is found by LHCb at 
about 2.8 CeV with = 3“ in the D^'k~ mass distribution, which is close to the mass of D 3 
state predicted in Ref. 
interaction. 


48|, 2800 ~ 2850 MeV. In the present work we investigate the AT-multi-p 


Taking FCA to the Faddeev equations, the KNN interaction studied in Refs. (50|,1^ has proved 
accuracy when dealing with bound states, and obtains consistent results with the full Faddeev 
equations evaluation without taking FCA [s^l, or a variational calculation with a nonrelativistic 
three-body potential model 
calculations in Refs. 



which is also conhrmed by the recent results with new Faddeev 
Even though there is a different results claimed in Ref. 56| on 


the KNN system, the work of Ref. ^ clarihed the different kinematical between them using 


two different approaches in their investigations, the Watson approach and the truncated Faddeev 
approach. A further study of the KNN system is done in the recent work of 58|, which has 


investigated the Kd scattering length with the hrst-order recoil correction using the non-relativistic 
effective field theory approach. A narrow quasibound state of 3500 MeV in the DNN system is 
predicted in Ref. [59( by both FCA to Faddeev equations calculation and the variational method 
approach with the effective one-channel Hamiltonian. Therefore, in the present work, we use the 
FCA to the Faddeev equations to investigate the A'-multi-p interaction. There are some possible 
K excited states with strangeness S' = ±1 and aligned with increasing spin number in the Particle 
Data Croup (PDC) [^, such as Ari(1270)(l’'') or Ari(1400)(l+) or Ari(1650)(l’''), Ar2(1580)(2“) 
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or i^2(1770)(2-) or K 2 ( 1820 )( 2 -) or i^2(2250)(2-), K3{2^20){3+), 7^4(2500)(4-). Some of these 
states still need more confirmation in the future experiments. This is the motivation of the present 
work, to understand their structure and properties theoretically. Applying the FCA to the Faddeev 
equations, there should be resonances or bound states in the two-body subsystem. For the iF-multi- 
p systems, the basic two-body subsystems are the pp and Kp interactions. Based on the local hidden 


gauge Lagrangians [6ll464l|. the two-body pp interaction is studied in Ref. [65!] with the coupled 


channel approach of ChUA and found that a pp quasibound state or molecule could be associated 
to the /2(1270) found in the PDG [^. With the on-shell Bethe-Salpeter equation of ChUA, and a 
chiral Lagrangian, the two-body Kp interaction is studied in Ref. [ 66 |] and dynamically produced 


the iFi(1270) resonance. Furthermore, the iFi(1270) state is reinvestigated with detail using ChUA 


to analyse the experimental data in the later work of Ref. 67]]. Thus, the resonances /2(1270) in 


the pp interaction and Ari(1270) in the Kp interaction are the needed clusters in the formalism of 
the present work. 

In the next section, we will first present the formalism of the FCA to the Faddeev equations. 
Then, in the following section, the resonances /2(1270) and Ari(1270) are reproduced dynamically 
with the ChUA in the pp interaction and the Kp interaction respectively. Our investigation results 
are shown in Section HVl Finally, we finish with conclusions. 


II. FORMALISM 


For the three-body interaction as Faddeev suggested in Ref. 
T-matrix can be written as a sum of three partitions, 

T = T^ + T‘^ + T^, 


35l ]. the scattering amplitude of 


( 1 ) 


where partition amplitude T* (z = 1, 2, 3) includes all the possible interactions contributing to 
the three-body T-matrix with the particle i being a spectator in the last interaction. But, if there 
are resonances (or bound states) as clusters appearing in the two-body subsystem interaction, for 
example the cluster coming from T^, thus, we can assume that a cluster is formed by the two 
particles (named as particle 1 , 2 ) and is not much modified by the interaction of a third particle 
(particle 3) with this cluster. Therefore, assuming the cluster as the fixed center of the three- 
body system, we can take the FCA 3a, |42 h 46I| to the Faddeev equations. Then, the T^ partition 
amplitude contributes to the cluster for the FCA, and the FCA multiple scattering of the third 
particle with the components of the cluster is taken into account. Thus, we can rewrite the Faddeev 
equations of Eq. ([T]) easily (which is also developed by the ChUA), 


Ti —ti + tiG'oT 2 , 

( 2 ) 

T 2 =t2+ t2GoTi, 

( 3 ) 

T = Ti + T2, 

( 4 ) 


where T is the total three-body scattering amplitude, as depicted in Fig. [H From this figure, we 
can see that the Faddeev equations under the FCA are first a pair of particles (1 and 2) forming 
a cluster, and then particle 3 interacts with the components of the cluster, undergoing all possible 
multiple scattering with those components. Thus, the two partition amplitudes Ti and T 2 sum 
all diagrams of the series of Fig. [1] which begin with the interaction of particle 3 with particle 1 
of the cluster (Ti), or with the particle 2 (T 2 ). Finally, the scattering amplitude T is the total 
three-body interaction amplitude that we look for. The amplitudes ti and t 2 represent the unitary 
scattering amplitudes with coupled channels for the interactions of particle 3 with particle 1 and 
2 , respectively, which should be taken into account the isospin structure of the subsystem and 
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FIG. 1: Diagrammatic representation of the FCA to Faddeev equations. 


discussed in details for different cases in Section nYi Besides, Go is the propagator of particle 3 
between the components of the two-body subsystem, given by 


Gois) = 


1 


■ d^q 

2M^J ^ 


FR{q) 


1 


(s) — g 2 — mg -h i e ’ 


(5) 


where q^{s), the energy carried by particle 3 in the rest frame of the three particle system, is given 
by 


q«(s) = 


s + mg — 


( 6 ) 


with m 3 the mass of the third particle and Mr the mass of the cluster, and FR{q) is the form 
factor of the cluster of particles 1 and 2. The form factor of the cluster should be consistently with 
the theory used to generate the cluster. This requires to evaluate the wave functions in the ChUA 
reproduced the cluster, which has been done in Refs. j68l-l7Cll| for 5-wave bound states, 5-wave 
resonant states and states with arbitrary angular momentum, respectively. Since in the present 
cases, the generated clusters are 5-wave bound states, we only need the expression of the form 
factors for the 5-wave bound states, given by (69|, 


^ ^ ^ /|p|<AMp-g|<A' ^ 2uJi{p) 2 u}2{p) Mr - UJi{p) - Ul2{p) 


1 


1 


1 


Af = 


2 ui{p- q) 2u}2{p- q) Mr - uji{p- q) - U2{p - q) ’ 

[ _ ^ _ ^ ^ ^ 

J\p\<A' ^2a;i(p) 2u;2{p) Mr - LOi{p) -U2{p)^ 


(7) 


( 8 ) 


where cuj = \Jq"^ + mf (i = 1, 2, m* the mass) are the energies of the particles 1,2. We use a cut 
off A' to regularize the integrals of Eqs. ([7]) and ([ 8 ]) (also for Eq. ([5]), discussed in Section HV]) , 
which is the same as the one used in the loop function of the two-body interaction to reproduce 
the cluster 47, Thus, no free parameters are involved. Note that, in the present work, we care 
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about the dynamics close to the threshold and thus the method of cutoff is acceptable 

Taking the normalization of the field theory 7^ which has different weight factors for the 
particle fields, we must take into account how these factors appear in the single scattering and 
double scattering and in the total amplitude [^. In all the present cases, the cluster (also 
particles 1 and 2) is a meson and the scattering particle (the third particle) too, which are only 
related to meson fields. Thus, following Ref. 7^ we write the S matrix of single scattering. Fig. 
□ (a), 

= — Ri(27r)^ 6{k + kn — k' — k'^ 

11111 


^2^^ = — R2(27r)^ 5{k + kR — k' — k'^ 


1 1 


(9) 


( 10 ) 


where, k, k' {kR, k'^ are the momentum of initial, final scattering particle {R for the cluster), 
Wi, are the energies of the initial, final particles, V is the volume of the box where the states are 
normalized to unity and the subscripts 1, 2 refer to scattering with particle 1 or 2 of the cluster. 
Next, the double scattering diagram. Fig. [I](b), is given by, 

111111 


5 ( 2 ) ^ _ i{2T:)^S{k + kR-k' -k'j^)^' 


d^q 

(27r)’ 


FR{q)- 


1 


Q 


02 




— g 2 _ Tie 

where FR{q) is the cluster form factor that we have discussed above, seen in Eq. (|7|l. 

Similarly, the full S matrix for scattering of particle 3 with the cluster can be written as, 

S = -iT{2TT)^6{k + kR-k’-k’R)x 11111 


( 11 ) 


122 ^2ids V^ioR 


( 12 ) 


Now, we can see that for the unitary amplitudes corresponding to single-scattering contribution, 
one must take into account the isospin structure of the cluster and write the ti and t 2 amplitudes in 
terms of the isospin amplitudes of the (3,1) and (3,2) systems. In view of the different normalization 
of these terms by comparing Eqs. Q, (11011 . (jlip and (|12ll . we can introduce suitable factors in the 
elementary amplitudes. 


2Mr 
— — - ti, 


2mi 


~ 2Mr 
t2 - - - t2, 

2772-2 


(13) 


where we have taken the approximations. 


But, when the cluster, the particles 1 and 


y/2uJi \/2mi ’ 

2 , one of them is a baryon, the factors in Eqs. dSl) and (fTdIl . 2Mr and 2mi should be replaced by 

1 for taking the baryonic field factor approximation 

obtained 


I 2Mb 

2Eb 


1. Finally, we sum all the diagrams. 


rp rp , rp h + t2 + 2 tl t2 Gq 

1 = 1^+12 = -^- r r ^2 - 


(14) 


^ We use a common cut-off method to regularize the propagator. Since it violates gauge invariance, it is a effective 
way to remove infinities from perturbative calculations, of which the problem is addressed in Ref. 0 within 
the effective Lagrangian. Furthermore, respecting Lorentz invariance, the cut-off regularization can recover the 
symmetry (references therein). 
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When ti = t 2 in some cases, it can be simplified as, 


T = 


2ti 


1 — tiGo 


(15) 


Note that, the FCA to Faddeev equations are particularly suited to study system with the subsys¬ 
tem bound or even loose bound, as discussion in Refs. [51 


5ll |. and have some limitation on the 


case when the cluster of the two particles in the subsystem is excited in the intermediate states 
(more discussions seen in Ref. [751] for the study of (pKK system). 

From the S matrix of single scattering, Eqs. ©, m, and the full S matrix, Eq. m , we 
should note that the arguments of the amplitudes Ti{s) and ti{si) are different, where s is the 
total invariant mass of the three-body system, and Si are the invariant masses in the two-body 
subsystems. The relationship between them is given by 


= m^ + mf + 


{Ml + mj - m2)(s - mi - M|) 


2M| 


, {i,j = 1,2, i / j). 


(16) 


III. TWO-BODY INTERACTION 

Using the Faddeev equations under the FCA, we first need bound states in the two-body sub¬ 
system as the cluster of the fixed center, and then let the third particle collide with the cluster 
and interact with the components of the forming cluster. As discussed in the introduction, for the 
subsystems of the A'-multi-p system, the two-body pp and pK interactions are studied in Refs. 

and [^ . We briefly summarize the works of them to reproduce the resonances /2(1270) and 
Ari(1270), and obtain the two-body scattering amplitudes of the subsystem (details given in the 
appendix). 


A. pp interaction 


In Ref. 65|], the pp interaction is studied with the local hidden gauge formalism [611464] and the 
ChUA with coupled channels, and dynamically produced the /2(1270) ( I{J^^) = 0(2'''''') ) state. 
From the local hidden gauge Lagrangians, the s-wave potentials of spin S = 2 for pp interaction 
are obtained in the sectors of isospin 7 = 0 and 1 = 2, 


rr“-"=^T.) = -4/-89y4|_i), (17) 

(18) 

where g = Mv/2fT^, with My the vector meson mass and /,r the pion decay constant. The 
scattering amplitude of pp interaction is calculated by the on-shell Bethe-Salpeter equation, 

= [l- V^G^]-W^, (19) 


where the kernel is a matrix of the interaction potentials and G^ a diagonal matrix of the loop 
functions (upper index I represents the specific isospin sector). 

Following the work of Ref. 65| (more details seen Appendix [^ , we also take into account the 
contribution of the box diagram with two pseudoscalar mesons in the intermediate state. We only 
consider the imaginary part of the box diagram contribution as the correction of the potential , 
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1.8e+06 



and neglect its real part which is very small. Note that we do not take these intermediate channels 
in the box diagram (more detail seen in Ref. 65]) as accounting for the coupled channels 0 , 0 - 


On the other hand, as done in Ref. 65||, we also consider the p mass distribution by replacing the 


loop function in the corresponding channel with its convoluted form. We obtain consistent results 
with Ref. 165(1, as shown in Fig. [21 where the structure of the resonance /2(1270) is found in the 


peak of the modulus squared of the amplitude. We have successfully reproduced the /2(1270) state 

(1=2 S=2) 

as the cluster of our procedure. Besides, the nonresonant amplitude tpp ’ ' is not shown in the 

figure, which is needed when we evaluate the two-body interaction amplitudes ti, t 2 , considering 
the isospin structure of the subsystem (discussed in Section iTVl) . 


B. pK interaction 


The pK interaction is investigated by the ChUA in Refs. [ 66 l . [67j. Following the work of 
Ref. [67(1 , we can dynamically generate the itri(1270) resonance in the interaction of pK and its 
couple channels, 4>K, uiK, K*r] and K*'k. From the local hidden gauge Lagrangians, the vector- 
pseudoscalar potential projected over s-wave can be written as 


Vij{si) — -gj2 ~ “I" “I" 

--{Mf-m]) , 


( 20 ) 


where and represent the masses of z (?) channel of the incoming (outgoing) particles, 

and the coefficients of Cij can be found in Ref. [67| (also given in Appendix|B]). Then, we can input 
these potentials into the the on-shell Bethe-Salpeter equation to evaluate the scattering amplitude. 


= [l + F^G^]-i(-F^)r-e', (21) 

2 

where is (1 + being a diagonal matrix {G^ as the normal loop function in Eq. (1191) i and 

e(e ') represents a polarization vector of the incoming (outgoing) vector-meson. As done in Ref. 
[67[ (details of the convolution of G^ given in Appendix [C|) , we also take into account the large 
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width of the vector mesons, and consider the convolution of the vector mesons as intermediate state 
in the loop function . In Fig. [3l we show our results for the modulus squared of which 

are consistent with Ref. 67|. From the peak position, we can see that the R'i(1270) is successfully 
reproduced in our work as the cluster in the FCA. In / = 3/2 sector, there are only two coupled 
channels, pK and iF*7r, which is no resonance appeared. But for considering the isospin structure 
of the subsystem, we also evaluate the amplitude which is not showed in the figure. 


IV. RESULTS 

In the present work, we study the iF-multi-p interactions. In the former section, we have 
reproduced the resonances /2(1270) and Rri(1270) in the pp and pK two-body interactions, which 
are the clusters of the FCA to Faddeev equations for the three-body interaction. Therefore, based 
on the possible clusters in the two-body interaction discussed above, the possible cases for the K- 
multi -/3 interactions are listed in Table [U and explained as follows. For the three-body interaction, 
we have two options: (i) particle 3 = K, cluster or resonance R = f 2 (particle 1 = p, 2 = p) and 
(ii) 3 = p, R = Ki {1 = p, 2 = K). For the four-body interaction, we can extrapolate the FCA 
ideas and also have two cases: (i) 3 = / 2 , R = iFi (1 = p, 2 = K) and (ii) 3 = Ki, R = f 2 
(1 = p, 2 = p). If we find a new resonance in the four-body interaction, assumed as K^, thus, 
there are also two cases for the five-body interaction, {i) 3 = K, R = (1 = / 2 , 2 = / 2 ) and (ii) 

3 = p, R = K‘i (1 = / 2 , 2 = Ki). For the six-body interaction, (i) 3 = ATi, R = /4 (1 = / 2 , 2 = / 2 ) 
and (ii) 3 = f 2 , R = K 3 (1 = f 2 , 2 = Ki). We show our investigation results for all these cases as 
below. 


A. Three-body interaction 

In the three-body interaction, there are two possible structures: K — / 2 (pp) and p — Ki{pK), 
which means (i) 3 = AT, R = /2 (1 = p, 2 = p) and (ii) 3 = p, R = R'l (1 = p, 2 = K). To 
evaluate these scattering amplitudes, we need as input the ti and t 2 amplitudes of the (3,1) and 
(3,2) subsystems, G = t 2 = tpK for K — / 2 (pp) and ti = tpp, t 2 = tpK for p — Ki{pK). We can 




TABLE I: The cases considered 
TO 3 -I- Mr (unit: MeV). 

in the R-multi-p 

interactions. 

Threshold 1: m 3 + mi + m2] threshold 

particles: 

3 

R (1,2) 

amplitudes 

threshold 1 

threshold 2 

Two-body 

P 

K 

tpK 

1271.0 

(1270.0)k, 

P 

P 

tpp 

1551.0 

(1275.1)y, 

Three-body 

K 

P 

/2 {PP) 

Ki ipK) 

Tx-h 

Tp-Ki 

2046.5 

1770.6 

2045.5 

Four-body 

Ki 

h 

/2 {PP) 

Ki (pK) 

Txi-h 

Tf2-Ki 

2821.0 

2546.1 

2545.1 

Five-body 

K 

U (/2/2) 

Tk-u 

3045.7 

2513.5 

P 

Rs (/ 2 R 1 ) 

Tp-Ks 

3320.6 

(2855.5) 

Six-body 

Ki 

f2 

/4 (/2/2) 

R 3 (/ 2 R 1 ) 

Tki-u 

Th-Ki 

3820.2 

3288 

(3355.5) 


calculate them in the two-body op and pK interactions, which are discussed in the former section, 
following the work of Refs. [67[. But, note that, in their work, the dimensional regularization 
scheme is used for the loop functions. To evaluate the form factor of the cluster, we need a cutoff A', 
which is the same as the one qmax used in the loop function for the two-body interaction, discussed 
in Section El As discussed in Ref. [78|], we can compare the value of the G function at threshold 


using the dimensional regularization formula 79l| w ith the one of the cut off which can be taken 


from Ref. 281] or the analytic expression in Ref. 80(|. Then, equivalently to the parameters in the 


dimensional expression, we obtain Qmax = 875 MeV for the /2(1270) cluster and Qmax = 1035 MeV 
for the Ail (1270) cluster. In fact, we do not introduce any free parameter. 

With the values of these two Qmax for the cutoff of A', then using Eqs. © and dH]), we can 
evaluate the form factors of the clusters, /2(1270) and Ari(1270), shown in Fig. 0](Left), where we 
can see that when q —)• 2A', Fji{q) —>• 0. Therefore, for the cutoff of the Gq function, seen in Eq. 

we choose as 2A' from the constraint of the form factor. Note that, the form factor of Eq. ([7|) 
is only valid for the cluster of bound state, which is the starting point of the EGA, since only the 
wave functions of the components of the bound state lead to the ordinary form factor but not for 
the cluster of resonance [691] (this is the limitation of the EGA, as discussed in the formalism). In 
fact, we also have two more form factors for the clusters, the and states, seen in Table [H 
Since the form factor often reduces fast above a certain momentum (q ~ 2A'), as shown in Fig. 
m (Left), we can safely choose the same cutoffs as the /2 and Ki for the ones of the /4 and 
respectively to avoid introducing any new free parameters, seen in Fig. |4] (Right). Furthermore, 
even though we change the value of the determined qmax a bit, for example 10%, if it is still in 
the nature value 33||, the shape of the form factor will not be changed, seen Fig. U thus, our final 
conclusion will not be changed, and only the strength of the amplitude will be a little differences. 

We have mentioned in the formalism that we should take into account the isospin structure of 
the subsystems for the two-body amplitudes ti and ^ 2 - For the first case of K — f 2 {pp), the cluster 
of /2 resonance has isospin 1 = 0. Therefore the two p mesons are in an / = 0 state, and we write 
it in terms of the physical basis components 


Iw >'“■“>= ^(l(l.-l)> +l(-l.l)> -l(0.0)>), 


( 22 ) 


where |(1,—1) > denote \{ll,l‘l) > which shows the Iz components of particles 1 and 2, and 
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FIG. 4: The form factors of /2(1270) and iFi(1270) (Left), /4(2050) and 7^3(2090) (right). 


\pp means \pp Then, the third particle is a meson taken |/| >= >, obtained 

Tl}:^fl = <K-pp\i\K-pp>^"2^"^^ 

=(< < pp\m) ^isi + 132 ) i\K >(^ 5 ) 0\pp >( 0 - 0 )) 

(4i + ^ 2 ) 

(1,-1) >+|(-l,l)>-1(0,0) 


< ^1 ® ^( < (i,-i)l+ < (-1,1)1- < (0,0 


1 

L'2> 


1 

■3 L 


3, 1 ,3 1, 


2 A l^ , 
3 ^ 2 ’ ~ 2 ’ ^ ~ 


2 ,31,, 

3 ^ ^2’2^’°l 


(23) 


/I A l^ , 
" Vs ^ ^ 2 ’ 2 ^’°' 


^31 


1.3 3, , /I,,3 1, ^ 2.A 1, , 

1 ^ 2 ’ 2 ^’" > +V 31 ^ 2 ’" 2 ^’ > +V 3l^2’"2^’ > 


/2,,3 1, /1,T 1, 

"V3l^2’2^’°^"V3l^2’2^’°^ 


1 

+ 31 


1 ,3 1, , /2 A lx , 

3 + ^2’“ 2 ^’^^ + V 3 + ^2’“ 2 ^’^^ 


3 3 

+ <(2>2)’"V- 


2 ,31,, I A lx , 

3 + ^2’2^’°l " Vs + ^2’2^’°l 


^32 


1,,3 1, 

3 2’ ~2 ’ ^ + 


+ 


2,,1 1, ,,3 3, /2,,3 1, /1,,1 1, 

3 ^^2’“ 2 ^’^ > + 1 ( 2 ’ 2 ^’“^ > “V 3^^2’2^’° + “V 3^^2’2^’° + 


where the notation of the states followed in the terms is |(|,|),—1 >= l(A^, Iz^)Az > for ^3i, and 
l(A^, Iz^)Az > for ^ 32 - Finally, we obtain the amplitudes combining with the isospin structure 

A = Ak = ^(2 t3r+" + 4r+"), A = A, (24) 

where pK scattering amplitudes with isospin I = 1/2, 3/2, and have been evaluated 

in Subsection IIIIBl 

But for the second case of p — Ki{pK), the isospin structure relationship is different. Now, the 
isospins of p and Ki are Ip = ^ and Ik^ = thus, the total isospin of the three-body system are 
two cases I total = IppK = ^ or Aotai = IppK = f • Therefore we have 

\pKi ><ii)= WK ><i' 4 )= / 11 ( 1 ,-i) > - 1 ( 0 . i) >, 

Ipifi \ppK ><t' 4 )= 1 ( 1 ,-i) > + yl 1 ( 0 , 1 ) >, 
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(25) 




















where we have taken the third isospin component Iz = \ for convenience. Then the \pK > states 
inside the Ki for the Iz = —h and Iz = +h are given by 




ipK >(«)=,/? 1(1.-i)> - 


i(o.h>. 


(26) 


For the two possibilities, using Eqs. (I2^ and (l26|) and performing a similar derivation of Eq. (I23I1 . 
we obtain 


p(^=l/2) . 
-p-Ki ■ 

p(/=3/2) . 

-p-Ki ■ 


2 .r_, 


^1 — tpp — 2 ^31 


1=0 


t2=tpK = ^(8 4 ^/^ 




); 


^1 — tpp ~ g ^31 ) t2 — tpK — g (5 ^32 ^ +4 t^2 

and 




(27) 


where pp interaction amplitudes with isospin I = 0) 2, t^i and are evaluated in Subsection 
nil A[ and the ones of pK interaction, ' and ' , are given in Subsection IIII Bl 

In Fig. Owe show our results of the modulus squared of the amplitude for There is a 

clear and sharp peak around 1770 MeV, which is close to the threshold of K — f 2 and similar to a 
cusp in the case of the p'KK system 81|. Therefore, this peak in would be affected by the 


threshold effect, and the width of this resonance structure has large uncertainties. There would be 
the cusp corresponded to a real resonance in some cases, like ao(980) (more discussions about 
the states appearing near the threshold can be found in the recent works of [s^ . 83|). In Fig. [6]we 
show the results of (left) and (right). From the 1444^ results, we find that 

there is a clear peak around the energy 1707 MeV with a width about 113 MeV which is about 
340 MeV below the p — Ki threshold and the ppK threshold. Because of the large width of the p 
meson, for a system with two p mesons, the large bindings in the present case will be acceptable. In 


PDG 


60|, the iir2(1770) of = 2“ strangeness state is the mass of 1773 ± 8 MeV and the width 
186 ± 14 MeV. But, from the analysis of the Kuj spectrum in the reaction K~p —)• K~ujp, the 
work of Ref. obtains its mass as 1710 ± 15 MeV and width 110 ± 50 MeV, which is consistent 
with our results. Our results are also consistent with the other experimental results 8^-87]. Thus, 
considering the uncertainties in our study (which will be discussed later), the peak appearing in the 
corresponds to the it'2(1770), which would be the p — Ki molecular state in our model. 

The strength of the peak of is about 25 times smaller than for |444^’ ^^^is, we could 

not expect a state structure in Fig. [5l even though the two-body interaction of pK is not strong as 


pp by comparing the results of Fig. [2]and Fig. [3l From the results of |T 


-/=3/2|2 


m 


Fig. [6] (right), we 


can see that there is a clear resonant structure about 2100 MeV with the strength 25 times smaller 
than that of \Tp~_ 44^ in the left figure, which is a little above the p — iFl(1270) threshold, and 
also show a dip in the threshold. We are looking for the lowest lying states bound in the RT-multi-p 
system, and hence, we could not expect a new state in the 1444^ results. 


^ Since the complicate situation in the definition of the second Riemann sheets in the multi-body interaction, we 
have the difficulty to extract the masses by the poles in the second Riemann sheets, which is different from the 
two-body interaction, more discussions seen in Ref. [^ . Thus, we extract the mass and the width from the 
scattering amplitude in the first Riemann sheet, shown in the figures. 
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FIG. 5: Modulus squared of the TK-f 2 - 



FIG. 6: Modulus squared of the Tp-Ki scattering amplitudes. Left: Itotai = Right: Itotai = §■ 


B. Four-body interaction 


For the four-body interaction as shown in Table HI we also have two possibilities: (i) particle 
3 = /2, cluster R = Ki [1 = p, 2 = K), or (ii) particle 3 = Ki, resonance i? = /2 (1 = p, 2 = p). 
Because the isospins of the two cluster are 7/2=0 and Iki = |, the total isospin of the four-body 
system is only Itotai = Using the FCA formalism as discussed in Section |TT] for the four-body 
interaction, for the first option, /2 interacting with the Ki, we need to evaluate the amplitudes 
ti = tf^p = which has been done in Ref. 411] and t 2 = tf^K = TK-f 2 ^ which has 

been calculated in the former Subsection IIV A1 Similarly, for the second case, Ki collides with the 
/2, and the amplitudes ti = t 2 = txip = Tp-Ki have been evaluated in the former Subsection 
IIV Al Note that the three-body amplitude Tp^Xi should be also written in terms of the isospin 
structure as discussed in the Section m and its amplitudes of isospin components evaluated in the 
last Subsection IIV Al Since the isospins of both the Ki and K are 7 = i, the isospin structure is 
similar to the case when the K collides with the f 2 - Thus, analogously to Eq. (|24p we have 


h = TpK, = -{2Tl^^/^ + Tl^^/^), 


h — ti- 


( 28 ) 


® where they used a similar formalism of Sec. m thus, using the formalism of Sec. El and the pp interaction 
amplitudes in Subsec. IIII Al we can reproduce their results. 
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FIG. 7: Modulus squared of the (left) and TKi-f 2 (right) scattering amplitudes. 


We show our results in Fig. [71 of which the left is and the right where the 

clusters are resonances Ki and /2 respectively. In the left of Fig. [71 there is a clear peak at an 
energy of about 2079 MeV, the width of which is about 249 MeV. We also hnd discontinuity at 
the position about 2550 MeV, which comes from the contribution of the /2 — Ki threshold effect. 
In the right of Fig. [71 we also find that there is a resonant peak around the energy 2091 MeV 
with a width of 230 MeV in the The strength of the peak of is at the same 


energy of the peak is a little less bound. In PDG, there 


magnitude as the one of 

is only one = S"*" iFa state found in the experiments, 7^3(2320), with a mass 2324 ± 24 MeV 
and width 150 ± 30 MeV. This 1^3(2320) state was just found in the reactions K^p —(Ap)p and 
K~p —?■ {Ap)p 89|. But, its mass is too far away from our results. Therefore, we hnd a new 
state with uncertainties, with a mass about 2079 ~ 2091 MeV and a width about 230 ~ 249 MeV. 


C. Five-body interaction 


As we expected before, we also find a new state in the four-body interaction in the former 
Subsection IIV Bl Thus, for the five-body interaction, there also are two options for the cluster, one 
of which is the /4 state found in the PDG and studied in Ref. [^, and the other one the resonance 
obtained in the four-body interaction above. Then following the idea of FCA and letting the 
third particle {K or p) collide with them, the two possibilities are (i) particle 3 = K, cluster ii = /4 
(1 = /2, 2 = /2), or (ii) 3 = p, R = {1 = f 2 , 2 = Ki). Since the isospin = 0 and Ik^ = 

for the first case, the total isospin of the K — f 4 system is only Itotai = but for the second option 
p — K^, the total isospin of this structure is I total = ^ or Itotai = Therefore, the situation of K 
interacting with /4 {K — /4) is similar to the three-body interaction discussed before, K colliding 
with /2 (7^-/2), and p — analogous to the one of p — Ki. Thus, in the first case, the K collides 
with the /4, and the amplitudes ti = t2 = tKf 2 = have been evaluated in Subsection IIV Al 

for the three-body interaction. For the second case, the p interacts with the which is similar 
io p — Ki in Subsection IIV Al thus, doing a similar derivation as done in Eq. (|23|) . we obtain 


1 / 2 ) . j. _ j. _ rpiI—1) , _ , _ rpl—ll‘2_ 

^ p-Ks ■ ~ ^Ph ~ -^31 ) *2 — tpKi — 132 5 

rp{I = 3/2) _ , _ , _ rp{I = l) , _ i _ rpI=3/2 

^p-Ks ■ ~ W 2 — ^3^ , C2 — — 432 ; 


(29) 
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FIG. 8: Modulus squared of the Tx-f^ (left) and Tp-x^ (right) scattering amplitudes. 


is the amplitude of Tp_f^, which is the same as calculated in the Subsection 


r^I = 3l2 


where the T, 

IIV Bl reproducing the results of Ref. [^, and the amplitudes T^Zx^ and T^Zx^ have also been 
evaluated in Subsection IIV Al 

We show our results for the two cases of the hve-body interaction in Fig. [ 8 l The results of 
is shown on the left of Fig. [ 8 l where we can see a resonant peak around the energy 
2505 MeV with a width of about 32 MeV. The right of Fig. [ 8 ] is the results of \T^pZxl^ and 

\TpZx‘^'^- observe that there are clear peaks for both of them. A resonant structure in \TpZx^\'^ 
is found at the energy 2382 MeV with the width about 409 MeV, which is about 120 MeV more 
bound than the one of and much larger width too. But the strength of \T^pZx^\'^ is just 

one half smaller than the one Thus, the more bound energy and lager width in \TpZx^\‘^ 

come from the stronger pp interaction in the components of A'a state, where we can see from the 
Fig. [2] and Fig. [3l In experiments, one = 4“ state was also found in Ref. [^, A4(2500), with 
the mass 2490 ± 20 MeV and width around 250 MeV, which is not well confirmed in PDG because 
of the lack of more experimental information. This 1^4(2500) state is dynamically generated as a 
molecular state oi K — f 4 ^ in our results of with a mass 2505 MeV and a width about 

32 MeV. The predicted width of our results is eight times smaller than the one reported. We 
should admit that there are also uncertainties in our results, seen from the results of on 

the right of Fig. [ 8 ]and discussed later. For {Tj^Zx^l'^ there is also resonant structure at the position 
about 2636 MeV, the width of which is about 171 MeV, which is not found for any I = 3/2 K 4 
state in the PDG. Thus, within uncertainties, we predict a new K 4 resonance of isospin I = 3/2, 
with a mass about 2636 MeV and a width about 171 MeV. 


D. Six-body interaction 


From Table m if we predict a K 3 state, seen in the Subsection IIV Bl analogously to the five-body 
interaction, there are also two options of the cluster for the six-body interaction, the particle 
found in PDG and studied in Ref. [411] and the resonance predicted above. Under the FCA, we 


let a third particle of the composite resonance {Ki or / 2 ) to collide with them, having (i) particle 
3 = Ki, cluster R = /4 (1 = / 2 , 2 = fZ)., or (ii) 3 = / 2 , R = R 3 (1 = / 2 , 2 = K\). Since the isospin 
of the particles, If^ = =0 and Ix^ = Ik^ = ^nd that the total isospin of the six-body 

system is only Itotai = Thus, for the first case, the Ki colliding with the {K — / 4 ), we need 
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FIG. 9: Modulus squared of the Tk^-j^ (left) and (right) scattering amplitudes. 


calculate the amplitudes ti = t 2 = tKif 2 — which have been evaluated in Subsection IIVBI 

For the second case, the /2 colliding with the K 3 , we evaluate the amplitudes ti = tf^f^ = 
by reproducing the results of Ref. ^ within our formalism (we have already reproduced Tp_f^ in 
Subsection IIV B() . and t 2 = tf^Xi = taking from the results of Subsection IIVBI 


We show our results for the six-body interaction in Fig. [9l The left of Fig. 

( 1 = 1 / 2 ) 


His 


Since Ki collides with the the amplitudes ti = t 2 = tKif 2 = , as shown in Subsection 


IIV B1 the amplitude for a test we also take ti = t 2 = tii'i/j = 


to 


evaluate the scattering amplitude again. By taking ti = t 2 = tii'i/j = we find a peak 

around the energy 2558 MeV with a large width of about 531 MeV. For the results of taking 
h = t 2 = tKif 2 — \ there also is a resonant peak in the position about 2670 MeV, the 

width of which is about 543 MeV. We can see that this test just gives us some uncertainties for 
our results. The right of Fig. [9] is , where we can observe that there is not a clear peak 

at the position 2681 MeV. This peak looks like the resonant structure of /2 — K 3 which is less 
stable comparing with the one of Ki — /4. Since there is no particle found in PDG, there 
could be a new resonance which is more uncertainty from our results, with a mass of about 
2558 ~ 2681 MeV and a large width about 531 ~ 543 MeV or more. 


E. Discussions 


Now we would like to make a discussion about the uncertainties of our results. As discussed in 
Section imi we have considered the width of the vector mesons in the evaluation of the two-body 
interaction amplitudes, tpp and tpx, by taking into account the convolution of the loop function 
. On the other hand, for the iF-multi-p system in the present work, the third particle in many 
cases, seen in Table. HI is a vector meson or a resonance, and has a large width, which should be 
taken into account. Thus, as done in Ref. [4^, we can roughly consider the width of the third 
particle in the propagator function Go, Eq. ([5]), by replacing the term in the denominator ie with 
Or, we can do the convolution of Go function for considering more exactly the contribution 
of the width of the third particle as done in Ref. 90|]. But, as shown in the results of Refs. 


90|, the final results for the position of the generated states are not altered after considering 


the contribution of the width of the third particle, to which we assign some uncertainties on the 
width for the generated results. Furthermore, there is also certain width for the clusters, for 
example, the resonances of /2(1270) and Ari(1270) in the three-body interaction, which have a big 
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FIG. 10: Considering the convolution of U in the K — pp scattering amplitudes. 


width. As done in Refs. 48], l90[], we also would take into account roughly the contribution of the 
width of the cluster by replacing the mass of the cluster Mr in Eqs. ([7]) and (l8|) with Mr — iTii/2. 
In fact, the contribution for the width of the cluster is just a small effect on the masses and widths 
of the generated results (seen more discussions in Refs. 48, [^). For checking this contribution, 
we also can consider the convolution for the since Eq. (I16p is dependent on the mass of the 
cluster, seen in Fig. [10] taking K — pp scattering for example, where we can see that the effect of 
the width of the cluster is small for the position of the peak and just makes the width of the peak 
a bit larger. Therefore, in our present work, we ignore all these effects in our investigation, just 
keeping in mind that there also some uncertainties from considering the contribution of the width 


of the third particle and the cluster. Besides, according to the discussions in Ref. [9]|, there may 


be some uncertain effect in the denominator term of ti t 2 Gq {ti Gq) in Eq. (fT4]l (Eq. (fT5]l l. 

We would like to make some discussions concerning the off-shell contributions in our formalism. 
For the two-body interactions, we follow the works of Refs. ^ for the pp and pK interactions, 
where the on-shell Bethe-Salpeter equation was used. A short prove of the on-shell Bethe-Salpeter 


equation is given in Refs. 28], l29|, and a different justification is given in Ref. 33]], where a sub 


tracted dispersion relation is used and one can obtain the results of the on-shell Bethe-Salpeter 
equation when one neglects the contribution of the left hand cut. Besides, the on-shell approxima¬ 
tion and full off-shell effects are discussed in detail in the work of 9^ for the vrvr scattering and in 
the works of Bonn group 9^95]. In the present work, we investigate the R-multi-p interactions 
with the EGA to Faddeev equations, which are different from the full Faddeev equations with ChUA 
as done in Ref. 39|. One of the intrinsic differences between them is that the EGA to Faddeev 
equations based on the two-body on-shell amplitudes as input while the full Faddeev equations 
using the full off-shell amplitudes. Yet, this should not be seen as a draw back of the FGA. Indeed, 
in the full Faddeev equations calculations of Ref. [s^, the full off-shell amplitudes can be separated 
into two parts: an on-shell part, and a remnant off-shell part which are unphysical. In full Faddeev 
calculations these unphysical parts will be implicitly canceled by ad hoc three-body forces which 
are included in the calculations, details seen in Refs. 9a, l97|. This means that using the on-shell 
amplitudes in the Faddeev equations is a more appropriate way to deal with the three-body system 

. The 


than using the full-off shell amplitudes if no additional three-body forces are introduced [ 
use of the on-shell amplitudes in the FGA gets a strong support from these findings. Indeed, since 
the on-shell two-body amplitudes are used as input in the FGA to Faddeev equations, the cancella¬ 
tion of the off-shell contribution can not be hold as in the full Faddeev equations calculations 39|. 
Up to one-loop, there are only two diagrams, seen (b) and (fj of Fig. [H while having six diagrams 
in the full Faddeev equations calculations. Following Ref. [39|, we divide the on-shell and off-shell 
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parts for the two-body amplitudes, written ti{si) = tf^{si) +t°/^{si), thus, using Eq. m we have 
for one-loop 



d^q 


Fniq)- 


g ,02 _ ^2 _ J ^2 ^ g 








(30) 


where we can see that the off-shell parts of the two-body amplitudes can not be canceled since 
that we do not have the interaction force from the contact terms. But, one also can be see that, 
because of the constraints of the form factor F^{q) (seen Fig. 0 ]), the contribution of the off-shell 
parts can not be infinity. The results of Ref. [99(1 (their Fig. 6) show that the position of the peak 
move to higher less then 3% when one take a qi dependent to ti{si). Thus, the contribution of the 
off-shell parts are small which can be included in our uncertainties. 


V. CONCLUSIONS 

In the present work, we investigate the many body interactions in R-multi-p systems, using the 
formalism of the fixed center approximation to the Faddeev equations. We start from the two-body 
interaction of pp and pK with the combination of dynamics of the local hidden gauge Lagrangian 
and the coupled channel effect, to reproduce the resonances of /2(1270) and ii'i(1270) as the clusters 
for our formalism. We summarize our results in Table HH where we obtain some bound states in the 
R-multi-p interactions, as the findings of the three-body wave functions in a finite volume in Ref. 
[^ . In the three-body Kpp system, we dynamically generate the R2(1770) state in our formalism, 
obtained a resonant peak in the modulus squared of the scattering amplitudes around the position 
in 1707 MeV with a width about 113 MeV and explained its structure as a p —Ri(1270) molecular 
state. Continuing with the three-body Faddeev equations with the fixed center approximation, we 
extrapolate our formalism to the four-body interaction. We observe a new state, with a mass 
about 2079 ~ 2091 MeV and a width about 230 ~ 249 MeV, which would be a Ki — /2 molecular 
resonance and not found in the PDG yet. For the five-body interaction, we successfully generate 
the R4(2500) state in our results of with a mass 2505 MeV and a width about 32 MeV, 

even though our theoretically predicted width is smaller than the one of unconfirmed experimental 
results of about 250 MeV. Thus, we also explain the structure of R'4(2500) as a molecular state of 
K — f^. Besides, we find a new K 4 resonance of isospin I = 3/2 in the p — interaction, with a 
mass about 2636 MeV and a width about 171 MeV. Finally, analogously we predict a new state 
in the six-body interactions of the Ki — and /2 — with a mass of about 2558 ~ 2681 MeV 
and a large width about 531 ~ 543 MeV or more, with more uncertainties. We hope that in future 
experiments, our predicted states of 7^3(2080), 7^5(2670) (isospin 7 = 1/2), and 774(2640) (isospin 
7 = 3/2) are found. 
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TABLE II: Our results of the A'-multi-p interactions (units: MeV). 


Interaction 

interactions 

Results (mass, width) 

Exp. Res 

Pridictions 


K- f 2 {pp) 

cusp 

- 

- 

Three-body 

p-Ki {pKY=^l'^ 

(1707, 113) 

A2(1770) 

- 


p-Ki {pKy-^/^ 

- 

- 

- 

Four-body 

h - Ki {pK) 

Ki - /2 {pp) 

(2079, 249) 

(2091, 230) 

— 

A 3 (2080) 


K-U ( 72 / 2 ) 

(2505, 32) 

A4(2500) 

- 

Five-body 

p-A 3 (/2Vi)"='A 

(2382, 409) 

- 


P-K 3 (/2Ri)"=3A 

(2636, 171) 

- 

A4(2640) 

Six-body 

Ki - fi ( 72 / 2 ) 

72 - R 3 (72 Vi) 

(2558-2670, 531-543) 
(2681, > 500) 

— 

As (2670) 


Matter (Hadron Physics 3, Grant No. 283286) under the Seventh Framework Programme of the 
European Union. 


Appendix A: Contribution of the box diagram and the convolution of the loop 


As mention in Ref. [65|, we should take into account the contribution of the box diagram to 
the potential . The main contribution is vrTr-box diagram for pp interaction, written 


T/(2’^d-0,S-2). X _ 

^box{pp) ys)-c>v , 


(Al) 


where is given by 


ISvr^ 

1 


f 


_ , ^ _ . /^ r^max „ „ 1 

={V2g)‘^ ' ,2 




I 


1 


(A2) 


fc? + ? — 2a; + ie A:? — ? — 2a; + ie — 2a; + 




le 


where to = ^ y/s = + k^- For the cut off in the integration, we took a natural 

size, q'max — 1200 MeV. Beside, F{q) is a form factor for an off-shell tt in each vertex, in the case 
of TTvr-box, taken 


m = 


h? — ml 

A2 _|_ ^2 ’ 


(A3) 


with A = 1300 MeV. 

The element of the loop function in the matrix in Eq. (11911 is only two meson loop, read 

d^q 


Giiis,mi,m2) = i 


1 


1 


(27r)^ q^ — nil + ie {P — qY — ni^+ie' 
which could be regularized by the cut off, obtained 

/•gmaa; ^2^1^ a;i-Fa;2 


Gjj(s,mi,m2) = 


'q (27r)2 a;ia;2[(P’°)^ - (wi UJ 2 Y + ie] ’ 


(A4) 


(A5) 


where a;* = yjq'^ + , the centre-of-mass energy {P^Y = •S; qmax stands for the cutoff, of 

which we took qmax = 875 MeV for the pp interaction. It is notice that we also could do the 
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TABLE III: Coefficients Cij in eq. in the I = 1I2 sector. 



(j)K 

ujK pK 

K*r] 

AT^tt 

(jjK 

0 

0 

0 

-i/i -Vi 

ojK 

0 

0 

0 

V3 

2 

v/3 

2 

pK 

0 

0 

-2 

3 

2 

1 

2 

K*r] 

[3 

V 2 

^/3 

2 

3 

2 

0 

0 

K*tt 

fs 

V 2 

^/3 

2 

1 

2 

0 

-2 


TABLE IV: Coefficients Cij in eq. (l20l) in the I = ij2 sector. 



pK K*tt 

pK 

1 

1 

K*tv 

1 

1 


calculation of Eq. (lASp with the analytic expressions in Refs. |3ll . [80|. But we need take into 
account the convolution due to the p mass distribution by replacing the G function as follows: 


Gpp{s) -^2 


-{mp+2rpf Y 

drrql — 


Im 


1 


ttJ m\ — rrfip + irhiT {rhi) 


r(mp+ 2 Vp)'^ 

hmp-2VpY 




Im 


ttJ rh2 — m'i + im2T{rh2) 


Gpp{s,rh‘l,ml), 


(A6) 


with 


N = 


f{rnp+2rpf 2 

/ dmi ( — — ) Im 

l{mp- 2 rp)^ 


ttJ + i?fiir(mi) ’ 


r(mi) = r, 


mf — 4m^N 3/2 


TUp — 




(A7) 

(A8) 


where Ep = 146.2 MeV and Gpp(s, , m^), is given by Eq. 


Appendix B: Coefficients of the potential 

The coefficients Cij in Eq. ()20p for the I = lj2 sector are given and tabulated in Table HIH 
In the I = 3/2 sector, the relevant channels are pK and K*Tr and the coefficients are given in 
Table HYl 


Appendix C: Convolution of loop function 

Taking the on-shell factorization, the loop function for two mesons is expressed as a function of 
the energy s, seen Eq. (|A4p . To remove the ultra violet divergence of the loop function, we follow 
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the dimensional regularization scheme 


Gii{s,Mi,mi) = 


1 


IOtt^ 
Ql 


, , , M? mf — M? + s, mf 




+ ^ [ln(s - {M^ -mf)+ 2qiy/s) 
V'S 

+ln(s + {Mf - mf) + 2qiyfs) 
-ln(-s + {Mf - mf) + 2qiy/s) 
-ln(-s - {Mf - mf) + 2ql^/^)] } , 

with a momentum qi determined at the center of mass frame 


m = 


^[s - {Ml - mz)2][s - {Ml + mz)2)] 

2Vi 


(Cl) 


(C2) 


where /i is a scale parameter in this scheme. The finite part of the loop function is stable against 
changes of /r due to the subtraction constant a{^) which absorbs the changes of fj,, where we take 
the following parameter set chosen to reproduce /Ci(1270) in Ref. 671 

/X = 900 MeV, a(/x) = —1.85, / = 115 MeV. (C3) 

But, considering a finite width of the vector mesons in the loop function, as done in Ref. (^ . 
the effect of the propagation of unstable particles is taken into account in terms of the Lehmann 
representation, which is done by the dispersion relation with its imaginary part 




ImD(sy) 
s — sy + ie’ 


(C4) 


where Sth stands for the square of the threshold energy. Now the spectral function is taken as 

1 


ImZl(sv) = Im 


sv — My + iMyTy 


(C5) 


where the width Ty is assumed to be a constant physical value. Substituting Eqs. ()C4p and 
into the original loop function Eq. m, we have 


Gii{s,Mi,mi) = ^ [ 

L; J(j 


{Mi+ 2 ri)^ 


(M,-2ro2 

-l|lm 

TT 


dsvGii{s, mi) 


1 


sy — Mf + iMiTi 

where Gu is given by Eq. m, and the normalization for the Ith component 

1 


|■( A ^^+ 2^^)2 / I 

Cl = dsv X I — 

J(Mi-2ViY V 


TT 


Im 


sy — Mf + iMiTi 


(C6) 


(C7) 


with mi, Ml, E;, the mass of the pseudoscalar meson, mass of the vector and width of the vector 
respectively. Replacing Gu by Gu in Eq. m, we include the width effect of vector mesons. In 
the present case, we only do the convolution for the p and K*. 
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